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ABSTRACT 


With the reduction of funding for sea ice reconnaissance flights, the 
National/Naval Ice Center needs to capitalize on the improvements in satellite 
technology. Imaging sensors such as AVHRR, DMSP/OLS, SSM/I and RADARSAT are 
used to detect the presence of sea ice, but with the exception of SSM/I and RADARSAT, 
clouds are a major obstacle to viewing the surface. With NASA’s development of the 
Moderate-resolution Imaging Spectroradiometer (MODIS) and MODIS Airborne 
Simulator (MAS), there is finally a sensor capable of using multi-spectral techniques to 
detect the presence of clouds. 

A group at the Space Science and Engineering Center (SSEC), University of 
Wisconsin — Madison lead by Dr. Steve Ackerman has developed a cloud mask for 
MAS/MODIS. The technique determines a level of confidence that a given pixel is clear 
based on a series of multi-spectral tests. By combining the confidence level from all 
tests, it is possible to detect the presence of clouds at different altitudes in the 
atmosphere. Threshold optimizations are described in this thesis for the Tg(11pm) and 
Tp(3.9um) - Tp(1 1pm) tests from Ackerman et al. (1997). In addition, the Tp(11um) - 
Tp(12um) test is removed. These modifications are based on daytime analysis of several 
MAS cases and a limited number of MODIS cases. 

Subjective analysis shows the modifications greatly improve the detection of 
clouds over cold polar oceans where sub-pixel ice may be present or water temperatures 
might falsely indicate clouds. The number of Cloudy pixels (<0.66 clear confidence 
level) for a given scene was increased 12.1% on average for MAS cases. The NPS cloud 
mask also classified two times more Probably Clear and Undecided pixels than the 


original mask due to greater sensitivity to thin, small clouds. 
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INTRODUCTION 


A. MOTIVATION 


With the reduction of funding for sea ice reconnaissance flights, manned polar 
observation stations, and polar in situ measurements, the National/Naval Ice Center 
(NIC) needs to capitalize on the latest advances in remote sensing. Clouds can be 
characterized by either their brightness temperature or reflectance characteristics. 
Previous cloud masks used simple combinations of visible and IR wavelengths to detect 
the presence of clouds. Cloud detection is difficult in polar ocean regions where sea ice 
may have the same brightness temperature and/or reflectance of the surrounding clouds. 
Detection becomes further complicated when clouds are thin along their edges and are 


almost indistinguishable from the sea ice. 


The Moderate Resolution Imaging Spectroradiometer (MODIS) and the MODIS 
Airborne Simulator (MAS) provide an opportunity for multi-spectral cloud detection 
approaches using many more wavelengths than previously available. Multi-spectral 
cloud detection will take advantage of the many differences in the spectral properties of 
sea ice and clouds. As a result, future automated sea ice analysis algorithms will more 
correctly identify sea ice properties with minimal cloud contamination. 


B. HISTORY 


Operational analysis of sea ice using remote sensing assets date back to the 
1970’s. Analysts manually analyzed sea ice in imagery with grease pencils and 
eventually transferred the analyses to paper charts. In the early 1990’s, on-screen sea ice 


analysis became possible due to increased computer power and products were 


disseminated via the World Wide Web. The future objective of operational sea ice 
analysis is an automated process with minimal manual intervention. Automated analysis 
requires a cloud mask that is sensitive to the unique conditions present in polar ocean 


regions. 


Satellite sensors are available that can penetrate clouds to derive surface 
information without the need of a cloud mask. One is a passive microwave sensor called 
Special Sensor Microwave Imager (SSM/I) on the Defense Meteorological Satellite 
Program (DMSP). SSM/I techniques usually use 85GHz imagery to identify sea ice and 
cloud liquid water content. SSM/I is not the perfect ice detection instrument. Microwave 
techniques currently provide low-resolution observations (ca. 25km) and have difficulty 
identifying sea ice concentrations when water pools on top of the sea ice. A second 
active technique uses the Synthetic Aperture Radar (SAR). For example RADARSAT 
uses a 5.3GHz radar to penetrate clouds and view the surface. SAR is an ideal tool for 
operational sea ice analysis but is restrictive due to the cost per scene to the U.S. 
Government. 


C. OBJECTIVE 


The objective of this thesis is to optimize a cloud mask for MODIS that is 
optimized for the polar ocean region and used by NIC and other polar scientists. The 
work described here is based on a cloud mask that was originally developed by the Space 
Science and Engineering Center at the University of Wisconsin - Madison. The 
optimized cloud mask will assist in the operational sea ice analysis by identifying clouds 
that might be mistaken as sea ice or their shadows mistaken as leads. The ultimate goal is 


the initialization of an automated sea ice analysis algorithm. 
2 


D. ORGANIZATION 


Chapter II outlines background information about the MAS/MODIS cloud mask. 
Chapter III describes the methodology for optimizing the MAS/MODIS Cloud mask. 
Chapter IV shows MAS examples of the NPS cloud mask adjustments. Finally, Chapter 
V shows a MODIS example of the NPS cloud mask and Chapter VI outlines the 


conclusions of this thesis and offers recommendations for further studies. 
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Il. BACKGROUND INFORMATION 


A. PLATFORMS AND SENSORS 


Two platform/sensor packages are employed by the cloud mask algorithm. First 
is the ER-2 and MAS research sensor, which is used to evaluate channel properties. 
Second is the Terra Satellite and MODIS operational space-borne sensor. 

Is ER-2 and MAS 

MAS is a modified Daedalus Wildfire scanning spectrometer that is flown on a 
NASA ER-2 aircraft. It has been used in many experiments including WINCE (January 
— February 1997) and FIRE/ACE (May — June 1998). It provides 50 band, 50m 
resolution from 20km flight altitude in 50 spectral bands range from 0.523um to 
14.383um wavelengths (Table 2-1). Band spectral ranges can change depending on the 
experiment. 

2: Terra Satellite and MODIS 

The Terra satellite, formally AM-1, is the flagship of the Earth Observing System 
(EOS). EOS is a major program under the National Aeronautics and Space 
Administration’s (NASA) Earth Science Enterprise (ESE) with a primary mission to 
observe the Earth’s surface. Terra was launched from the Western Test Range, 
Vandenburg Air Force Base, California on December 18, 1999 at 13:57 EST (Herring, 
2000). It is in a 10:30am local equator crossing, descending, Sun synchronous orbit in 
attempt to maximize the view of the surface while minimizing cloud cover. The orbit 
will be adjusted to ensure a complete surface coverage every 16 days with a life span of 6 


years (NASA, 1999a). 


Table 2.1. MAS Instrument/Platform Specifications 


Aircraft at 20 km altitude 
View ; 
Pixel Spatial Instantaneous 
Resolution Field of view . 


Pixel per Scan Line | 716 (roll corrected) | Roll Correction 


Scan Rate 6.25 scans/second Data Channels 50 channels via digital 
recorded 


Spectral Bands 


Bandwidth Bandwidth 
0.523 — 0.569 2.875 — 3.039 


























0.627 — 0.68 3.031 — 3.195 
0.679 - 0.721 3.202 — 3.358 
0.72 - 0.762 3.338 — 3.507 
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One of the five instruments aboard TERRA is MODIS. It has 36 bands ranging 
from 0.4um to 14.4141m wavelength. The bands are designed to observe land and ocean 
temperatures, land surface coverage, aerosols, water vapor, temperature profiles, and 
fires. Table 2-2 lists the various spectral bands, their associated spatial resolutions and 
primary use. 


B. MAS/MODIS ARCTIC CLOUD MASK 


The MAS/MODIS Cloud Mask is described in the algorithm theoretical basic 
document (ATBD) (Ackerman et al., 1997). The physical rationale and detail 
specifications of the cloud mask are documented in the ATBD. This thesis will only 


describe the important factors that effect cloud detection in polar ocean regions. 


Figure 2.1 shows MAS and MODIS bands overlaid on the transmittance 
properties for ozone, water vapor and carbon dioxide for a mid-latitude summer 
atmosphere calculated by MODTRAN 3.5 (Acharya et al., 1993). In general, cloud and 
surface observations are made in spectral regions with minimal atmospheric absorption. 
However, absorption bands can be used to limit the radiance from lower altitudes. This 
can help identify upper-level clouds. 

1. General overview of the Cloud Mask 

It is evident that the MAS and MODIS sensors have much in common in the 
exploitation of transmission and absorption features. In some cases they used different 
wavelengths to achieve the same result. As an example, MAS uses a 1.884m band to 


exploit low and mid level strong water vapor absorption to highlight the presence of high, 


Table 2.2. MODIS Characteristics (NASA, 1999b 


Bandwidth Spatial Resolution 
adir 


Tand/Cloud 
Boundaries 
[and/Cloud 
Properties 
6 | 108-1652 um | 500m 
[8 _0.405-0.420um | 100m 
(9) 0438-048 um | «000m 




































Ocean Color/ 
Phytoplankton/ 
Biogeochemistry 




















Atmospheric 
Water Vapor 













Surface/Cloud 
Temperature 


























Temperature 
Water Vapor |_27_| 6535-6895 um | 1000m 
pee Ea ister [tom 
Temperature [32 [| 11.770 12.270um [100m 
Cloud Top 
Altitude 

P36 [14.085 —14.385 um 


OplxOIp uoqies 
pue 1odea JoyeMm *9Uu0ZO JO} Somodoid sdue}}IWIsUe OY) oe plejIoAGQ ‘oloydsoune ay} ynoqe uOneULIOJUT UTe}gGO 
0) sotiedoid o1ioydsoune snouea oy) Wojdx9 syuoWIN.sUI osoy) s1aYyM diysuoyejes Aouonbajy - spueq SIGOW/SVW 7 ound 


(win) yybuajanepA 


spueg SIGOW —°— 
spueg SVN —e— 


Jequny pueg 
soueyiwsuel | 








thin cirrus clouds. To accomplish the same task, MODIS uses the water vapor absorption 
band at 1.38u1m. 

There are five major groups of multi spectral tests in the mask. Each group of 
tests are designed to detect common cloud features ranging from low clouds to high 
cirrus clouds. The number of tests within a group and the number of groups used in a 
particular case vary depending on the time of day, the ecosystem, presence of snow or sea 
ice, amount of perceptible water, and several others. Once an individual pixel’s 
characteristics have been described (e.g. daytime, ocean, polar, ice), tests within various 
groups will be preformed along with specific thresholds for that pixel’s characteristics. 
All tests within a group and all groups have the same weight as any other test or group 
used for a pixel. 

a. Calculation of Confidence 

Ackerman et al. (1997) developed a simple linear scaling of thresholds. 
This technique allows for a degree of uncertainty to be included in each test. Figure 2.2 
depicts the linear function used by the cloud mask. There are three points; alpha (q) - 
any value below this point is considered cloudy, beta (B) - this is the 50% clear 
confidence value (the pixel can be either clear or cloudy) and gamma (jy) - any value 
above this point is considered clear. Each test within a group is given a value between 0 
and 1 depending on its threshold. A value of 0 is an indication of a high probability of 
cloud, able a value of | is an indication of a high probability of a clear condition. The 


lowest confidence value from the tests within a group becomes the value for the over all 


group. 
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The overall clear confidence value for the pixel, Q, is computed for the 


pixel as the N" root the product of group confidence values: 


Q= uTT G, , (2-1) 


where G; is the value for a group and N is the number of groups used in the test. This 
calculation is conservative in the estimate of a clear sky condition and will err on the 
cloudy side for any given pixel. For example, if a pixel were in a daytime, ocean, polar, 
sea ice location, the cloud mask would use four groups for the test. Confidence values 
for the four individual groups might be 0.82, 0.66, 1, and 1. Therefore the over all 
confidence would be [(0.82)(0.66)(1)(1)]* = 0.86 for that pixel. Table 2.3 shows Q 
values for the four confidence levels. The pixel in this example would be assigned as 


Undecided classification. 


CLEAR CONFIDENCE 


CLOUDY 


TEST VALUE 





oO 
Figure 2.2. Clear confidence level linear function used by the cloud mask. 
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Table 2.3. Confidence Values 


 — Gonfitence | 0h 


Undecided Q> 0.66 — 


Cloudy or Obstructed | Q < 0.66 





a Inputs 
The cloud mask requires some information about the scene before it can 
determine if a pixel is clear or not. These inputs determine the number of groups and the 
number of tests that will be preformed within each group. Some of the variables that are 
important for the polar ocean portion of the mask are day/night, ecosystem, and sea ice 
detection. 
a. Day/Night 
The day/night calculation is based on the solar zenith sun angle (8). When 
6 > 85° it is considered day time (Ackerman et al., 1997). 
b. Ecosystem 
The ecosystem is determined by the pixel’s location within various index 
files and may be manually selected. They are: 


li Global USGS 1 km resolution land/sea tag 


os Loveland | km resolution N. American ecosystem 
a 10 minute resolution global Olson ecosystem 
4. Global Loveland 1 km resolution Olson ecosystem 


For the MAS cases, the Global Loveland 1 km resolution Olson ecosystem 
file was used to identify the ecosystem. The ecosystem is identified as either being 


desert, land, coast or ocean. 
iW 





Cc. Sea Ice Detection 

There are challenges to overcome when it comes to the detection of sea 
ice. First, sea ice detection is difficult when there are low clouds in the scene due to their 
similar radiance and thermal signatures. Second, relatively thin sea ice detection is 
difficult because polar ocean waters might have similar radiance values. 

The MODIS cloud mask has two ways to detect the presence of sea ice. 
First, Special Sensor Microwave Imager (SSM/I) analysis at 25km resolution, available 
from either the National Center for Environmental Predictions (NCEP) or the National 
Snow and Ice Data Center (NSIDC), is used to identify the presence of sea ice. The 
inputs from NCEP and NSIDC are given equal weight in the mask for the determination 
of sea ice extent. Unfortunilty, SSM/I has diffuculty in the summer when liquid water 
pools on top of sea ice and might incorrectly indicate an open water condition. 

Second, an internal means for idenfitying sea ice was developed by Hall et 
al., (1998) (herafter referred to the Hall Test) and implimented in the mask in an attempt 
to determine the pixel-scale extent in sea ice. In the Hall Test, when the brightness 
temperature at 1luum, Tg(11Ium), is less than 280°K, the surface is assumed to be cold 
enough to have snow or sea ice present. The difference between snow or sea ice depends 
on whether the ecosystem is either land/coastal or ocean. Next, the Normalized 
Difference Snow Index (NDSI) test detects sea ice when its thickness is greater than 
10cm. This works because sea ice absorbs 1.6211m wavelength radiance. Therefore, 
when sea ice is present, visible reflectance, Possum, Will be greater than 1.62)4m 


reflectance, Pi.62ym. In addition, sea ice has a high reflectance at 0.88um wavelength 


(Po.sgum) and Possum can be used to identify its presence. Hall et al. (1998) defines the 


NDSI as 


NDST = Pessim ~ Pisa (222) 
Possum + P6240 
When sea ice is present, Tp(1 1pzm) will be below 280°K and the NDSI will be computed. 
If the NDSI > 0.4 and Possum > 0.1 then sea ice is possibly present. 

Finally, a test is needed to ensure that the assumed sea ice is not a high, 
cold cloud. Two bands are included in the Hall Test that are sensitive to high clouds. 
The atmosphere has a low tranmittance at 134m due to CO, absorption and at 1.38um 
due to water vapor absorption. In both cases, radiance from the lower troposphere is 
removed by absorption and only radiance emitted or reflected from the upper troposphere 
contributes to the radiance measured by the sensor. If Ta(13pm) < 220°K and Pi.38um > 
0.035, the Hall test assumes that high cloud is present in the pixel and is responsible for 
the false-positive assignment of sea ice. 

When applying the Hall Test, it is important to remember that polar ocean 
water has a very low albedo compared to that of the highly reflective arctic sea ice. 
Complications arise when dealing with “New Ice” (Frazil, Grease, Slush or Shuga) that 
has a thickness between 0 and 10 cm. There might not be enough contrast in reflectance 
(Possum) to adequately separate the sea ice from water. In addition, as sea ice melts, 
water will pool on top of the floe and lower the albedo so the pixel will appear to be 


water when, in fact, it is sea ice. 
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For the operational MODIS cloud mask, sea ice is present only when 
either the NCEP or the NSIDC sea ice product and when the Hall Test detects the of sea 
ice. MAS only uses the Hall Test to detect the presence of sea ice. 

3. Groups and Test 

The five groups within the cloud mask are comprised of various tests. The 
number of groups and tests are determined by the input values. The groups that are 
important to the ocean polar regions are described below. 

a. Group I 

The objective of Group 1 is to identify thick high clouds based on infrared 
Tg characteristics. This is accomplished in the ocean polar regions with Tg(11um) and 
Tg(13.911m) wavelengths tests. The group has a potential difficulty with low clouds that 
have small thermal contrast between the cloud and the surface (Ackerman et al., 1997). 

(1) Tg(11p1m) Test. This is a simple IR cloud test that can obtain 
information about the surface or cloud top temperature. It is also important when dealing 
with the detection of snow covered sea ice as described above. In the open ocean when 
Tg(11um) channel is expected to be less than 270°K ,the mask assumes the pixel is too 
cold to be open ocean and determines that it has failed the clear sky condition. In the 
polar ocean regions, where there might be subpixel sea ice present, there is a need to 
lower the Tg(1lpum) threshold values. The magnitude of this modification will be 
discussed in the next chapter. In general, simple thresholds should vary with location, 
season, and time of day. 

(2) Tg(13pum) Test. This test relies on the absorption 


characteristics of CO2. Clouds being sensed by a CO, band depend upon the atmospheric 
is 


attenuation in that band and the altitude of the cloud. The spectral band at 13.9um 
provides good sensitivity to the relatively cold regions of the atmosphere. Only clouds 
above 500 hPa will have strong contributions to the radiance at the top of the atmosphere 
observed at 13.91um. Lower regions of the atmosphere and the earth’s surface should 
provide negligible contributions to the radiance at 13.9 um wavelength (Ackerman et al., 
1297): 

b. Group 2 

The objective of Group 2 is to detect thin clouds missed by Group 1 using 
infrared Tg differences. Tp difference techniques using MODIS will be very sensitive to 
thin clouds if the surface emittance and temperature are adequately characterized. Given 
a surface at 300°K and a cloud at 220°K, a cloud emittance of 0.01 increases the infrared 
window Tg by 0.5°K. Since the anticipated noise equivalent temperature for many 
MODIS infrared window spectral bands is better than 0.05°K, the cloud mask should be 
very sensitive to thin cloud (Ackerman et al., 1997). 

(1) Tri-Spectral Technique. Ackerman et al. (1997), suggests the 
use of a tri-spectral technique with Tg(8um) - Tp(11um) and Tp(1 1pm) - Tg(12um) tests 
for thin cirrus clouds when precipitable water (pw) is greater than 1.0cm. A polar 
atmosphere is relatively dry and pw is less than 1.0cm. Therefore, only the Tg(11pm) - 
Tp(12um) test is used to detect thin cirrus clouds. As will be discussed in the next 
chapter, MAS data suggests that this test might not be viable in daytime polar ocean 


regions. 
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(2) Tp(3.9um) - Tg(11lum). This test is used to identify fog and 
low clouds. Clouds at 114m wavelength, independent of water droplet size, are very 
close to being black bodies, with emittance approximatley equal to 1. Clouds at 3.9um 
wavelength are not black bodies and emittance is a function of water droplet size. Fora 
normal water droplet size of 10um, the emittance is approximately 0.85. 

During the day, the Tp(11p1m) is composed of only energy that has 
been emmitted from the cloud and is very close to the actual cloud top temperature. The 
Ta(3.9um) during the day is composed of both energy that is emmitted and reflected from 
the cloud. This will raise the sensed temperature above the actual cloud top temperature. 
Therefore, a positive difference in the Tg(3.9um) - Tg(11um) test is an indication of 
liquid water clouds during the day. 

At night, negative values for Tg(3.9um) - Tp(1lum) indicate a 
liquid water cloud. Tp(3.9um) is only a function of emmitted energy and will be less 
than the actual temperature since the emittance is less than one. Therefore, small 
Tp(3.9um) - Tp(1 lum) differences are observed when clouds are optically thick due to 
reduced emmitance at Tg(3.9um). 

c. Group 3 

The objective of Group 3 is to detect low clouds by using solar reflectance 
and reflectance ratio tests and is especially sensitive to thick, low clouds. It complements 
Group 1 because it is sensitive to thick, low level clouds while Group 1 has difficulty 


with low clouds that have small thermal contrast between cloud and background. This 


1: 


technique works well when there is a high contrast in the reflectance between the surface 
and the cloud, for example, clouds over water. 

(1) Possum. This is a single band threshold test that does well for 
discriminating bright clouds over dark surfaces (e.g., stratus over ocean) and poorly for 
clouds over bright surfaces (e.g., sea ice). The reflectance test is adjusted for viewing 
angle and the test is eliminated when a sun glint condition exists for the pixel (Ackerman 
et al., 1997). 

(2) Ratio of Po.g7um/Po.c6um. Many earth surface features are less 
reflective at visible wavelengths than they are at near-IR wavelengths. Clouds do not 
exhibit these same variations in reflectance to any great degree (Ackerman et al., 1997). 

d. Group 4 

The objective of Group 4 is to detect high thin cirrus clouds during the 
day. Due to the fact that most of the atmosphere’s water vapor is located in the low to 
mid troposphere, most of the reflectance from the surface, low, and mid level clouds, 
should be attenuated prior to reaching any space borne sensor when sufficient water 
vapor (0.4 cm precipitable water) is in the atmosphere. Gao et al. (1993) developed a 
technique to detect the presence of thin cirrus clouds by a strong water vapor absorption 
band in the near-IR at 1.38um (MAS - 1.88um) wavelength. This band receives large 
amounts of scattered solar radiance when thin cirrus clouds are present. Ackerman et al. 
(1997) placed this technique in to its own group because they felt “there was no test that 


has been found to be as sensitive to the presence of thin cirrus.” 
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é. Group 5 

The objective of Group 5 is similar to those of Group 2 in that it is 
designed to detect thin clouds. However, in this case the focus is on thin cirrus clouds at 
night. In the polar ocean regions, a Tp(3.911m) - Tp(11pm) difference is employed to 
detect cirrus clouds and is relatively insensitive to the amount of water vapor in the 
atmosphere (Hutchison and Hardy, 1995). This test is greater than zero in 
semitransparent cirrus as sub-pixel warm features dominate the shortwave window 


radiances within a pixel. 
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Il. MAS/MODIS CLOUD MASK 


MAS offers the opportunity to test the validity of the cloud mask prior to the 
MODIS launch. Subjective analysis of the cloud properties are compared with cloud 
mask results. Modifications of tests and thresholds are then suggested to improve the 
cloud mask performance. 


A. OPTIMIZING THE CLOUD MASK 


MAS Flight #97-73 is a daytime flight from the FIRST ISSCP Regional 
Experiment - Arctic Cloud Experiment (FIRE/ACE) flight of June 2 - 3, 1998 and is used 
as a case for cloud mask optimization in the polar ocean region. Figure 3.1 shows the 


over-ground track of the ER-2 aircraft. 
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Figure 3.1. | The ER-2 over-ground track is plotted for MAS Flight #98-73, June 2 - 3, 
1998. (Gumley et al., 2000) The two regions of the flight used for cloud 

mask optimization are identified as blue squares. 
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1. Original Cloud Mask 


Figure 3.2 is an annotated RGB image (Red: Band 20, Green: Band 10, Blue: 
Band 2) of track 17, center time of 01:04:29 UTC with Pt Barrow to the top of the image 
(South). The scene contains an area of mid to low clouds over a pack of mixed multi- 
year and first-year sea ice that dominates most of the scene. The Marginal Ice Zone 
(MIZ) is the transition area between pack ice and lesser concentrations of ice. Take note 
of the belting of sea ice jetting into water from the main pack ice. An area of 40 - 60% 
concentration of multi-year and first-year sea ice is located in the upper left hand corner 
of the scene. 

The original cloud mask classifies the scene as one of two types, Ocean Polar 
Snow Day or Ocean Day. This classification is very dependent on the Hall Test 
described in Chapter II. Table 3.1 lists the original group, test and thresholds for the two 
scene types. 


Table 3.1. Original Groups, Tests, and Thresholds 
Polar Ocean Snow Da Ocean Da 
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Figure 3.2. MAS Flight #97-73, Track 17, Center Time of 01:04:29 UTC is displayed 
as an RGB image where Red is band 20, Green is band 10, Blue is band 2. 
Clouds, multi-year and first-year sea ice and partial ice concentration area 
are noted. 


Figure 3.3 shows the original cloud mask correctly classifying most of the clouds 
as Cloudy. It also correctly classifies the clear region over the pack ice as Confident 
Clear. The cloud mask encounters difficulty along the MIZ, where it classifies the clear 
areas as either Undecided or Cloudy. It also has trouble in the water areas where the 
pixel is clear of clouds and sea ice, but the mask still classifies it as either Probably 
Clear, Undecided or Cloudy. 


a. Hall Test 


The Hall Test is very important in determining which scene type is used in the 


cloud mask. Since the entire scene is below 280°K for Tg(11um), the Hall Test is 


BIS 





initiated everywhere in the scene. Figure 3.4 shows the NDSI and Po xgsum tests for this 
scene. With the exception of a few cloudy pixels, the NDSI test classifies most of the 
scene as possibly being sea ice. The Posgum test checks the scene reflectance and 
classifies the polar water area as free of sea ice if the reflectance is less than 0.10 (green 


areas in Figure 3.4). 


) Confident Clear 
(Probably Clear 
Undecided 
C1Cloudy 





Freure 3.3. eR ate Glow Mask is. aoeted to MAS Flight 497. fe) lrack 17, Center 
Time of 01:04:29 UTC. Mask classifications are described in the legend. 


Figure 3.5 shows that the Hall Test correctly identifies the pack ice and 
much of the cloudy region. The test identifies some of the clouds as sea ice but as will be 
seen latter, are eventually classified as Cloudy when the final confidence value is 


computed. Finally, the floating sea ice in the partial ice concentration area is correctly 
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Figure 3.4. | The NDSI test (a) and the Possum test (b) is applied to MAS Flight #97-73, 
Track 17, Center Time of 01:04:29 UTC. The figures are colored 
according to the thresholds set for each test. The NDSI test classifies 
pixels with values greater than 0.4 as Probably Ice. The Possum classifies 
pixels with values greater than 0.1 as Probably Ice. 
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identified as sea ice. Therefore, the Hall Test properly classifies pixels and initiates the 


groups and tests that need to be conducted for cloud detection. 
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Fieure 3.5. MAS Flight #97-73, Track 17, Center Time of 01:04:29 UTC, Hall Test. 
Mask classifications are described in the legend. 


b. Individual Pixels 

Figure 3.3 reveals the greatest misidentification in this scene that occurred 
along the MIZ and the partial ice concentration area. Two representative pixels are 
selected to investigate this error. They are | and 2 pixels away, respectfully, from the 
main pack ice. Table 3.2 lists test results for a pixel that is classified as sea ice by the 
Hall Test and uses the groups and tests that correspond with Polar Ocean Snow Day. 
Table 3.3 lists test results for a pixel that is immediately adjacent to the first pixel and is 
classified as water and therefore uses the groups and tests that correspond with Ocean 


Day. 
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Table 3.2. Polar Ocean Snow Day Pixel Values for MAS Flight #97-73, Track 17, 
Center Time of 01:04:29 UTC, Scanline124304, Sample 305 
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Wable:3.3: Ocean Day Pixel Values for MAS Flight #97-73, Track 17, Center Time 
of 01:04:29 UTC, Scanline124304, Sample 306 
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When the cloud mask classifies the pixel as sea ice (Table 3.2), an overall 
confidence value of 1 is generated for the pixel. Even though the wrong scene type was 
initially selected, the cloud mask properly identifies the pixel as Confident Clear. When 
the cloud mask classifies the pixel as water (Table 3.3), the overall confidence value of 
0.36 is generated for the pixel. The correct scene type is selected but the cloud mask 
incorrectly identifies the pixel as Cloudy. The confidence values range from Possibly 
Clear to Cloudy when the surface is water and the confidence value is Confident Clear 


when the surface is sea ice. 
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The Hall Test may identify the pixel as water when a pixel possesses sub- 
pixel sea ice. The Po.ggum test threshold for the Ocean Day scene type is sensitive enough 
to detect the increased reflectance from the pixel due to subpixel sea ice. If the 
reflectance is too great for water, the clear confidence level from the Poggum test 
decreases. This occurs in the MIZ where sub-pixel sea ice occurrence is high. Tg(11i1m) 
and Tp(11lum) - Tg(12pm) tests also have low clear confidence levels (0.63 and 0.77 
respectfully) but are not as dependent on the presence of sub-pixel sea ice and are 
therefore better candidates for the threshold adjustments necessary for optimized cloud 
detection in polar ocean regions. 

c. Tzp(1 1m) 

The Tp(11p1m) test produced a 0.63 clear confidence level that the pixel 
was clear for the Ocean Day scene type case. Figure 3.6(a) shows the RGB image of the 
scene with the Tp(1 1pm) test result overlaid where the test was conducted. The clouds 
are not cold enough to be classified as Cloudy according to the original threshold, while 
the partial ice concentration area is not warm enough to be classified as Confident Clear. 

Surface temperature controls the value of Tp(11p1m) when corrected for 
surface emittance and given low water vapor content. Hall et al. (1998) used data from 
MAS images collected in April 1995, near St. Lawrence Island, Alaska, and estimated the 
mean seawater temperature to be 271.6°K using seawater emittance of 0.985. Based on 
direct measurement, the freezing point of polar ocean water was found to be 271.4°K. 
Therefore, any pixel can be classified as possibly sea ice if its surface temperature is less 


than or equal to 271.4°K, the freezing point of polar seawater. 
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Figure 3.6(b) is similar to 3.6(a) except the NPS modified threshold for 
Tp(11pm) test is applied. The threshold was lowered by 1°K from the original threshold. 
Polar ocean water is now warm enough to be identified as clear (white), while clouds are 
still cold enough to be identified as clouds since they are colder that the B-threshold (B = 
2727). 

When the NPS threshold is applied to the Tg(11tm) test and the cloud 
mask is run, the clear confidence level of the test pixel was increased to 0.79 from the 
previous value of 0.63. Figure 3.7 shows the original cloud mask (a) compared with the 
cloud mask with the modified NPS Tg(11pm) threshold (b). The clear confidence level 
increases for many of the pixels in the partial ice concentration area. The clouds are still 
identified as Cloudy. However, an unacceptable number of pixels are still classified as 
Undecided. 

d. Tz(11 pm) - Tg] 2pm) 

Figure 3.8 shows the Tp(11um) - Tg(12um) test colonzed for an average 
view angle and Tp(11lum) threshold. It is difficult to obtain any information about the 
scene from this test. Therefore, this test is removed from the Ocean Day scene type tests 
in the optimized cloud mask. 

Figure 3.9 shows the original cloud mask (a) compared with the cloud 
mask without the Tg(11um) - Tp(12um) test. There is significant improvement of the 
cloud mask in the partial ice concentration area showing fewer pixels as Cloudy or 
Undecided. Most of the water region has changed to Probably Clear with much of the 


noise is removed 
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Figure 3.6. 
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Original Tg(1 1pm) test threshold (a) and the NPS Tg(1 lpm) test threshold 
(b) are applied to MAS Flight #97-73, Track 17, Center Time of 01:04:29 
UTC and overlaid on a RGB images. The figures are colored according to 
the thresholds set for each test. The original threshold classifies values 
greater than 270°K as Probably Clear. The NPS threshold classifies 


values greater than 269°K as Probably Clear. 
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Figure 3.7. Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
Flight #97-73, Track 17, Center Time of 01:04:29 MAS Flight #97-73, 
Track 17, Center Time of 01:04:29 UTC. Mask classifications are 
described in the legend. 
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Figure 3.8. The Ta(11um) - Ta(12um) test is applied to MAS Flight #97-73, Track 
17, Center Time of 01:04:29 UTC. The scene is colorized for an average 
threshold. 
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Figure 3.9. | Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
Flight #97-73, Track 17, Center Time of 01:04:29 MAS Flight #97-73, 
Track 17, Center Time of 01:04:29 UTC. Mask classifications are 
described in the legend. 
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é. T2(3.9m) - Tg(1 1m) 

Once the new threshold was defined for Tg(11um) and the Tp(11um) - 
Tp(12pm) test was removed, other scenes indicated the Tp(3.94m) - Tp(11pm) test 
threshold might not be sensitive enough to detect low clouds. To be consistent 
throughout the algorithm, the Tg(11pzm) - Tg(3.9um) test in the Ocean Day scene type 
was changed to Tg(3.9m) - Tp(11um). The only difference in the test is a sign change. 

Figure 3.10 shows Track 9 of FIRE/ACE Flight #97-73, with a center time 
of 23:00:22 UTC and is an example where the entire scene is classified as sea ice on the 
surface and is dominated by thick low to mid-level clouds above. Dark shadows and 
semi-transparent clouds near the top of the scene might affect the performance of the 
cloud mask. Finally, along the bottom, multi-year sea ice is visible through the thin, 
broken low-level clouds. 

Figure 3.11 shows the original cloud mask. The original cloud mask over- 
characterized the scene as Confident Clear, especially in the dark shadow and semi- 
transparent area where high clouds are casting a shadow over low clouds. The original 
cloud mask also had difficulty detecting the thin, broken, low-level clouds. 

The cloud mask modified with the new Tg(1lpm) threshold and with 
Tg(1lum) - Tp(12p2m) removed was run but there were no changes in the detection of 
clouds. This is due to the fact that the entire scene was classified as sea ice and only used 
Polar Snow Day as the scene type. The Pigg and Tp(13.9m) thresholds correctly 


determined the scene contained no high clouds. The Tg(3.9um) - Tp(11pm) test is 
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therefore the only candidate for threshold modification in order to correctly classify the 


low and mid-level clouds apparent in Figure 3.10. 
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~~ i ~~ Thin, Broken Low-Level Clouds \ 
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Figure 3.10. MAS FIRE/ACE Flight #97-73, Track 9, Center Time 23:00:22 UTC, is 
displayed as an RGB image where Red is band 20, Green is band 10, Blue 
is band 2. Thick low to mid-level clouds dominate the scene with dark 
shadows and semi-transparent clouds near the top. Finally, along the 
bottom, multi-year sea ice is visible through the thin, broken low-level 
clouds. 


The Tg(3.9p1m) - Tg(1 1pm) test has a flipped threshold. The test indicates 
the scene is Confident Clear until it reaches the lower threshold value, and the confidence 
of the pixel being clear decreases as the test reaches the upper threshold value. At this 


point, there is zero confidence that the pixel is clear. 
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Figure 3.12(a) shows the orginal Tg(3.9um) - Tg(11um) test threshold. 
The original threshold was unable to identify many of the low clouds in the scene. It also 
failed to identify the clouds in the dark shadow and semi-transparent area. 
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Figure 3-11. Original Cloud Mask is applied to MAS FIRE/ACE Flight #97-73, Track 
9, Center Time 23:00:22 UTC. Mask classifications are described in the 
legend. 


Figure 3.12(b) is similar to 3.12(a) except the NPS modified threshold for 
T,(3.9um) - Tp(1 1pm) test is applied. The new threshold reduced @ from 11 to 5.25, B 
from 7 to 3.5 and y from 7 to 1.75. The new threshold improved the detection of the low 
level clouds in the thin broken and dark shadow, semi-transparent areas. The NPS 
Tp(3.9um) - Tg(11pum) threshold still had trouble where the dark shadows reduced the 
difference enough to be falsely considered clear. 

Figure 3.13 shows the cloud mask with the NPS modified Tg(3.9um) - 


Tg(1 1pm) test threshold added. The cloud mask detects more low and middle level 
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Probably Cloudy 


"NPS T,(3 Sum) -T,(14um) Threshold 


Figure 3.12. Original Tg(3.9um) - Tg(1 1pm) test threshold (a) and the NPS T,(3.9um) 
- Tp(1 1pm) test threshold (b) are applied to MAS FIRE/ACE Flight #97- 
73, Track 9, Center Time 23:00:22 UTC. The figures are colored 
according to the thresholds set for each test. The original threshold 
classifies values greater than 9.0 as Probably Cloudy. The NPS threshold 
classifies values greater than 3.5 as Probably Cloudy. 
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Figure 3.13. Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 

Flight #97-73, rack 9, Center Time “of 25:00:22 UTC) = wiask 
classifications are described in the legend. 
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clouds. Only the darkest shadow, semi-transparent area pixels are still classified as 
Confident Clear. 
B. NPS MAS/MODIS CLOUD MASK 

The observed weaknesses in the original cloud mask were revealed by two new 
scene type subroutines (Polar Ice Day and Polar Ocean Day). These scene types were 
developed after the removal of the Tp(1lum) - Tp(12um) test and by identifying the 
limitations in the thresholds for Tp(11um), and Tp(3.9um) - Tp(1 1pm) tests. These new 
scene types only run when the pixel is located in a daytime polar ocean region or when 
sea ice has been detected during the day. Table 3.4 lists the groups, tests and thresholds 


for these new scene type subroutines. 


Table 3.4. New Scene Type Subroutines. 
Polar Ice Da Polar Ocean Da 


Threshold Threshold 
(a, B, ¥) (a, B, 
266, 269, 272 
i T9(13.9um) 219,220,221 1 Bt LiL 
T,(13.9um) 219,220,221 
ag T,(3.9um)—Ta(lium) | 5.25, 3.5 1.75. ———— Ta(3.9p1m) — Ta(1 12m) 10, 8, 6 
8,7, 6.5 
0.95, 1.10, 1.15 
Po. | Baeitecm| Po.s6umm 0.95, 0.90, 0.85 
PS) wee ae 35, 


Figure 3.14 shows the final NPS cloud mask for Track 17, center time of 01:04:29 

































UTC with the new scene type subroutines. The clouds are better defined and cover a 
larger area. A few pixels over the cloud-free pack ice are now misclassified as 
Undecided due to the change in the Tp(3.911m) - Tp(11p1m) test threshold and appears to 


be related to noisy pixels. 


au) 


The new scene type subroutines greatly improved the classification of the partial 
ice concentration area. Most pixels are classified as either Confident Clear or Probably 


Clear - a direct result of the Tp(1lum) threshold change and the removal of the 


Tp(1 1pm) - Tg(12tm) test. 


The NPS cloud mask still has problems with the Po.ggum test when pixel coverage 
is Polar Ocean Day and sub-pixel sea ice is present. Under these conditions, the 
reflectance value rises above the threshold, but changing this threshold degrades the 
ability of the Polar Ocean Day scene type to correctly classify the pixel as Confident 
Clear when no sub-pixel sea ice is present. Misclassification is most evident along the 
MIZ where there is a large concentration of sub-pixel sea ice. The artifact can be used to 
identify the MIZ or other sea ice edge features. Unfortunately, the problems associated 
with sub-pixel sea ice contamination will increase as the resolution of the pixel decreases 
from 50m (for MAS) to 1km (for MODIS). On the other hand, much of the noise in the 
MAS IR bands will decrease because the comparable MODIS IR bands are expected to 


have better radiometric sensitivity and more accurate calibration. 
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Figure 3.14. Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
Flight #97-73, Track 17, Center Time of 01:04:29 UTC. Mask 
classifications are described in the legend. 
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IV. MAS DATA 


With the new thresholds in place, other MAS scenes were reviewed to test the 
optimization of the cloud mask. MAS data from the Winter Cloud Experiment 
(WINCE), First ISSCP Regional Experiment - Arctic Cloud Experiment (FIRE/ACE) and 
Winter Experiment (WINTEX) was used to verify the two new scene types. Subjective 
cloud analysis was done to verify the effectiveness of the new NPS MAS/MODIS cloud 
mask (NPS cloud mask). 

A. MAS FIRE/ACE FLIGHT #98-68, TRACK 6, CENTER TIME 21:04:36 


UTC 


The MAS FIRE/ACE Flight #98-68 (May 26-27, 1998), track 6, center time 


21:04:36 UTC is shown in Figure 4.1. 


0 4° [24:04:36 UTC | 
a 





FLIGHT 98-058 26-27 MAY 1998 a/c 809 FIRE-ACE (AMPR / MIR / CLS / NASBO / NISR / SSFR / HIS) 


Figure 4.1. © MAS FIRE/ACE Flight #98-68, May 26-27, 1998, Flight Track. (Gumley 
2000) 
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This case is an example of a relatively clear arctic condition with a high solar 
zenith angle (57.9°). Figure 4.2 shows the scene is dominated by sea ice with thin clouds 
along the right hand side. The Hall Test classified the entire scene as sea ice, therefore 
only the Polar Ice Day scene type was used to detect cloud. The Tp(13um) test detected 


no high clouds in the scene. 





Figure 4.2. FIRE/ACE Flight #98-68, Track 6, Center Time 21:04:36, is 
displayed as an RGB image where Red is band 20, Green is band 10, Blue 
is band 2. Sea ice dominated scene with a thin cloud along the right hand 
side. 

Figure 4.3 shows the original cloud mask. It classifies most of the scene as clear 


with a cloud along the right side of the scene. The striping in the scene is directly related 


to the noise in the MAS IR bands. 
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Figure 4.3. Original Cloud Mask is applied to MAS FIRE/ACE Flight #98-68, Track 
6, Center Time 21:04:36 UTC. Mask classifications are described in the 
legend. 


Figure 4.4(a) shows the original Tg(3.9um) - Tp(1lum) test threshold. The 
original threshold fails to detect a large portion of the thin cloud. Figure 4.4(b) shows the 
result of the NPS Tg(3.9um) - Tp(11pum) test. The NPS threshold detects much more of 
the thin cloud and possibly a very thin cloud along the bottom. The NPS threshold 
results in significantly more noise in the clear area, where it has lowered the final clear 
confidence value in many of the clear pixels. Again, the lowered values are an artifact of 


the MAS IR calibration. 


Figure 4.5 shows the original and NPS cloud masks. The NPS mask still 
classifies most of the scene as Confident Clear. The NPS threshold results in 
significantly Undecided pixels where it has lowered the final clear confidence value in 


many of the clear pixels. Identification of these Undecided pixels are an important result 
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Figure 4.4. 








Probably Cloudy 


Onginal T,(3.9,1m) - T.,(110m) Threshold 


Probably Cloudy 


NPS T,(3.9um) - T,,(11j1m) Threshold 


Original Tp(3.9.m) - Tp(1 1m) test threshold (a) and the NPS T,(3.9um) 
- Tg(11pm) test threshold (b) are applied to MAS FIRE/ACE Flight #98- 
68, Track 6, Center Time 21:04:36 UTC. The figures are colored 
according to the thresholds set for each test. The original threshold fails to 
detect a large portion of the thin cloud. The NPS threshold detects much 
more of the thin cloud and possibly a very thin cloud along the bottom. 
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Figure 4.5. 
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NPS Cloud Mask 
Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
FIRE/ACE Flight #98-68, Track 6, Center Time 21:04:36 UTC. Mask 
classifications are described in the legend. The NPS cloud mask classifies 
most of the scene as Confident Clear, and improves the detection of the 
thin cloud. 


47 





of the variable threshold techniques used in the MAS/MODIS cloud mask. Undecided 
pixels in Figure 4.5 identify regions of very thin clouds that are also apparent in Figure 
4.2. 

B: MAS FIRE/ACE FLIGHT #98-68, TRACK 16, CENTER TIME 23:31:38 


UTC 


Figure 4.6 shows the FIRE/ACE Flight 98-68, track 16, center time 23:31:38 
UTC. This is an example of a scene that contains both sea ice and high clouds. Clouds 
are casting a shadow over the sea ice that, at lower spatial resolution, could be 
misidentified as a lead. As the main cloud thins along its edges, it becomes increasing 


difficult to distinguish between cloud and sea ice. 


Figure 4.7 shows the original cloud mask. It identifies the main cloud mass but 
has difficulty along the edge. There is no indication of the thin clouds in the center of the 
scene suggested in Figure 4.6. The Hall Test classifies the entire scene as sea ice, 
therefore only the Polar Ice Day coverage subroutine is run. Again, Tg(13.9p1m) did not 


detect any clouds in the scene. 


Figure 4.8(a) shows the original Tg(3.9um) - Tg(1lum) test threshold. The 
original threshold detects much of the main cloud mass but subjective analysis suggests 
that this test is the reason why the original cloud mask is having difficulties with the thin 
cloud edges. Figure 4.8(b) shows the NPS T,(3.9um) - Tp(11um) test threshold. The 
NPS threshold is detecting the thin clouds that were missed by the original threshold. It 
has also detected two smaller clouds along the bottom right of the scene. Although this is 


a daytime case and these two clouds are in shadow, the emittance difference between 
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Figure 4.6. 


Figure 4.7. 


Shadow 





MAS FIRE/ACE Flight #98-68, Track 16, Center Time 23:31:38 UTC, is 
displayed as an RGB image where Red is band 20, Green is band 10, Blue 
is band 2. The scene contains both sea ice and high clouds, with clouds 
casting a shadow over the sea ice. 







|iConfident Clear 
‘DProbably Clear 
‘undecided 
‘Cloudy 


Original Cloud Mask is applied to MAS FIRE/ACE Flight #98-68, Track 
16, Center Time 23:31:38 UTC, Mask classifications are described in the 
legend. 
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Tp(3.9um) and Tg(1lpm) still reveals the presence of the clouds. Again, noise is 


introduced into the test by the MAS IR bands. 


Figure 4.9 shows the NPS cloud mask for this case. It detects the thin clouds 
along the edge of main cloud mass due directly to the threshold modifications in the 
Tp(3.9um) - Ta(1 1pm) test. It also detects the two smaller clouds along the bottom right. 
Noise in the cloud mask is generated by the Tg(3.9um) - Tp(1 11m) test. 
ay MAS FIRE/ACE FLIGHT #98-68, TRACK 4, CENTER TIME 20:11:05 


UTC 


MAS FIRE/ACE Flight #98-68, Track 4, Center Time 20:11:05 UTC, is a 
difficult case because sea ice is visible through semi-transparent clouds. The case is 
shown in Figure 4.10 where an optically thick cloud is casting a shadow on the clouds 
below it. Gravity waves are detectable in the thin, low-level clouds. Tg(13.94m) 


detected very little cloud and was not a factor in the ultimate confidence determination. 


Figure 4.11 shows the original cloud mask. It classifies much of the scene as 
Cloudy or Undecided. It also classifies a portion of the shadow and the lower right 


corner of the scene as Confident Clear. 
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Figure 4.8. 
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Shadow 


Clear 


NPS T,(3.9um) - T,(14p1m) Threshold 
Original Tg(3.9um) - Tg(1 1pm) test threshold (a) and the NPS Tg(3.9pm) 
- Tp(1 1pm) test threshold (b) are applied to MAS FIRE/ACE Flight #98- 
68, Track 16, Center Time 23:31:38 UTC. The figures are colored 
according to the thresholds set for each test. The original threshold did not 
detect the thin clouds, while the NPS threshold improved the thin cloud 
detection. 
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Figure 4.9. 
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NPS Cloud Mask 


Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
FIRE/ACE Flight #98-68, Track 16, Center Time 23:31:38 UTC. Mask 
classifications are described in the legend. (b) NPS Cloud Mask. The 


NPS cloud mask improved the detection of thin clouds. 
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Figure 4.10. 


Semi-Transparent Cloud 


Gravity Waves 
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MAS FIRE/ACE Flight #98-68, Track 4, Center Time 20:11:05 UTC, is 
displayed as an RGB image where Red is band 20, Green is band 10, Blue 
is band 2. Optically thick cloud is casting a shadow on the clouds below 
it. Sea ice is visible through semi-transparent clouds and gravity waves are 
detectable in the thin, low-level clouds 
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Figure 4.11. Original Cloud Mask is applied to MAS FIRE/ACE Flight #98-68, Track 
4, Center Time 20:11:05 UTC. Mask classifications are described in the 
legend. Much of the scene as classified as Cloudy or Undecided. A 
portion of the shadow as classified Confident Clear. The lower right 
corner of the scene, the cloud mask identifies the area as Confident Clear. 


Figure 4.12 shows where the Hall Test identifies sea ice. The pack ice in the 
scene and much of the multiyear ice floating in the open water is classified as sea ice. 


The test did not classify much of the optically thick cloud as sea ice. 


Figure 4.13 shows the Tp(3.9um) - Tp(11um) difference. Semi-transparent low- 
level clouds and gravity waves are resolved. The test also resolves the shadow from the 


thick liquid water cloud. 


When the Hall Test classifies the pixel as sea ice, the Polar Ice Day coverage is 
used. Figure 4.14(a) depicts the original Tg(3.9um) — Tp(1 1pm) test threshold. The test 


classifies the shadow area and much of the low level clouds as clear. The original test 
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resolves portions of the gravity waves in the low level clouds but the differences do not 


exceed the threshold necessary to classify the entire region as cloudy. 


| 
{Sea Ice 





Figure 4.12. MAS FIRE/ACE Flight #98-68, Track 4, Center Time 20:11:05 UTC, Hall 
Test overlaid on a RGB image. The pack ice in the scene and much of the 
multiyear ice floating in the open water is classified as sea ice. The test 
did not classify much of the optically thick cloud as sea ice. 


Figure 4.14(b) shows the NPS Tp(3.9um) — Tp(1 1pm) test threshold. The test 
now correctly classifies the shadow area as cloud. Much of the low level clouds that 
were missed by the original threshold are now identified as cloud. Finally, the gravity 


wave features are fully resolved by the test. 
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Figure 4.13. TT (3.9m) - Tp(11um) difference applied to MAS FIRE/ACE Flight #98- 
68, Track 4, Center Time 20:11:05 UTC. The color scale on the right is 
the range of temperatures found in this scene. Semi-transparent low-level 
clouds and gravity waves are resolved. The test also resolves the shadow 
from the thick liquid water cloud. 


Another Polar Ocean Day test is Tg(11pm). Figure 4.15(a) shows the original 
Tp(1 1pm) test threshold. The test correctly classifies the scene as cloudy. It has no 
trouble in the shadow area classifying it as cloud. Figure 4.15(b) shows the NPS 
Tg(1 lum) test threshold. It raises the clear confidence level where thin clouds are 


present over the sea ice but not to a level where it might be classified as Confident Clear. 


Figure 4.16 shows the original and NPS cloud mask. The NPS cloud mask 
correctly classifies the shadow due to the threshold modifications in the Tp(11um) and 
Tp(3.9um) - Tg(11um) tests. The low-level clouds are also correctly classified as Cloudy 


due to the modification in the Tg(3.9um) - Tg(1 1pm) test. 
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Figure 4.14. 
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“NPS Te(3.9um) - (114m) Threshold 
Original Tg(3.911m) - Tp(11"m) test threshold (a) and the NPS T,(3.9um) 
- Tg(1 1pm) test threshold (b) are applied to MAS FIRE/ACE Flight #98- 
68, Track 4, Center Time 20:11:05 UTC and overlaid on Tg(3.9m) - 
Tp(1 1pm) difference. The figures are colored according to the thresholds 
set for each test. The original threshold classifies the shadow area and 
much of the low-level clouds as clear. The NPS threshold correctly 
classifies the shadow and low-level clouds as clouds 
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lteure 4.15. 








Semi-Transparent Cloué 


Semi-Transparent Cla ic 


NPS T,(14j1m) Threshold 
Original Tp(11um) test threshold (a) and the NPS (11um) test threshold 
(b) are applied to MAS FIRE/ACE Flight #98-68, Track 4, Center Time 
20:11:05 UTC and overlaid on Tp(llum). The figures are colored 
according to the thresholds set for each test.. Both tests classify the scene 
as generally cloudy. 
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Figure 4.16. 
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Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 
FIRE/ACE Flight #98-68, Track 4, Center Time 20:11:05 UTC. Mask 
classifications are described in the legend. The NPS cloud mask correctly 
classifies the shadow and low-level clouds as Cloudy. 
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D. MAS WINCE FLIGHT #97-45, TRACK 6, CENTER TIME 20:18:05 UTC 


Figure 4.17 shows the ground track for Flight #97-45 conducted on February 6, 
1997. This flight in not in a polar ocean region but has lake ice and is in the area of 


responsibility of the NIC. 





FLIGHT 97-045 S FEBRUARY 1997 A/C 708 ASSO 
LAMBERT CONFORMAL PROJECTION: GP1 = 40-7 SP2 = 63.8 CH = -84.4 ROTATED BY 6.6 


17:21:14 70 22:42:38 UT SCALE 1:9.S53E+66 TIME TICK EVERY {6.06 MINUTES 


Figure 4.17. MAS WINCE Flight #97-45, Flight Track. (Gumley, 2000) 


Figure 4.18 shows an RGB image of MAS WINCE Flight #97-45, track 6, with a 
center time of 20:18:05 UTC. It is in James Bay between the coast and Akimiski Island. 
There is fast ice along the shoreline with a large floe in the center of the scene. Leads are 
located between the fast ice and the floe where no sea ice is apparent. Small clouds are 


casting elongated shadows on the sea ice due to the high solar zenith angle of 75.2°. 


The Hall Test classified the entire scene as either sea ice or snow with the 


exception of small areas that correspond to visible clouds. This included the leads that 
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appear to be open water in the visible and the IR bands. Figure 4.19(a) shows the NDSI 
test. With the exception of a few small areas where clouds are present, the NDSI value of 


the entire scene is greater than 0.4 and is therefore possibly sea ice or snow over land. 
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Figure 4.18 MAS WINCE Flight #97-45, Track 6, Center Time 20:18:05 UTC is 
displayed as an RGB image where Red is band 20, Green is band 10, Blue 
is band 2. Fast ice located along the shoreline with a large floe in the 
center of the scene. Leads are located between the fast ice and the floe 
where no sea ice is apparent. Small clouds are casting elongated shadows 
on the sea ice 


Figure 4.19(b) shows Po.sgum for this scene. The entire scene may contain either 
snow or sea ice because no Po.ggpm Value was less than 0.1. In the leads, where po ggum 1s 
expected to be less than 0.1, there were no values less than 0.12. This may be due to a 
thin sheet of ice in the leads that is thin enough to be transparent at visible wavelengths 


but increases Po.ssym enough to pass the test for sea ice. 
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Figure 4.20 shows the original cloud mask. It correctly identifies much of the 
scene as Confident Clear. The clouds were detected, but a subjective analysis suggests 
that the cloud mask could do a better job in depicting the size and shape of the clouds. 
With the lead classified as sea ice, the cloud mask correctly identifies it as Confident 


Clear. 


Figure 4.21 shows the Tp(3.9um) - Tp(11um) difference. This reveals the low 
level clouds in the scene. Figure 4.22(a) shows the original threshold as applied to the 
Tp(3.9um) - Tp(11pm) test. The original threshold detects some of the clouds. Figure 
4.22(b) shows the NPS threshold applied to the Tp(3.94m) - Tg(1 1pm) test . Subjective 
analysis reveals that the NPS Tp(3.9um) - Tg(11pm) test threshold increases the area of 
cloud detection. The clouds are better defined and better match the shadows in Figure 
4.18. When fast ice and sea ice are present, there is a decrease in the clear confidence 
level for this test. In leads, the test is negative, indicating that Tp(11um) is greater than 
Tp(3.9um), thus increasing the clear confidence level. Again, the MAS IR channels are 
noisy. 

Figure 4.23 shows the original and the NPS cloud mask. The new threshold for 
the Tg(3.9um) - Tp(11um) test has increased the number of pixels classified as Cloudy. 
Many pixels are classified as Undecided due to the noise in the Tg(3.9um) - Tp(11um) 
test, but these pixels still generate a confidence level of greater then 0.85 clear confidence 


level. 
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Figure 4.19 
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The NDSI test (a) and the Possum test (b) is applied to MAS WINCE Flight 
#97-45, Track 6, Center Time 20:18:05 UTC. The figures are colored 
according to the thresholds set for each test. The NDSI test classifies 
pixels with values greater than 0.4 as Probably Ice. The Possum classifies 
pixels with values greater than 0.1 as Probably Ice. 
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Figure 4.20 Original Cloud Mask is applied to MAS WINCE Flight #97-45, Track 6, 
Center Time 20:18:05 UTC. Mask classifications are described in the 
legend. 
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Figure 4.21 Tp(3.9pm) - Tg(11um) difference applied to MAS WINCE Flight #97-45, 
Track 6, Center Time 20:18:05 UTC. The color scale on the right is the 
range of temperatures found in this scene. 
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E. STATISTICAL ANALYSIS 


Table 4.1 shows the results of analyzing 18 tracks of data from three MAS flights. 
The average length of these tracks was 350 km and covered and average area of about 
13,000 km?. The NPS cloud mask classified 12.1% more pixels as Cloudy than the 
original mask. Subjective analysis of these tracks confirms that the new Cloudy pixels 
were indeed consistent with cloud features and not due to erroneous classification. 
Confident Clear pixels decreased by 35.2% in the NPS cloud mask indicating far fewer 
pixels are incorrectly identified as cloud-free. The NPS cloud mask generated about 
twice as many Probably Clear and Undecided pixels as the original cloud mask due to 
greater sensitivity to small, thin cloud features. 


Statistical summary of original and NPS cloud mask pixel classification. 


Original NPS Ghance 
Cloud Mask Cloud Mask 8 
Total Pixels 92388744 92388744 i 22 


% Confident Clear | 39781316 (43.0%) | 25769650 (27.9%) | -14011666 (-35.2%) 
% Probably Clear 1276156 (1.4%) 2444210 (2.6%) 1168054 (91.5% 
) 


Table 4.1 












% Undecided 4039600 (4.4%) 9200997 (10.0% 5161397 (127.8%) 
% Cloud 47291672 (51.2%) | 54973887 (59.5%) 5512242 (12.1%) 
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Figure 4.22 
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NPS T,,(3.9um) - T,(14m) Threshold 
Original Tp(3.9um) - Tg(1 1pm) test threshold (a) and the NPS Tp(3.9um) 
- Tp(1 lpm) test threshold (b) are applied to MAS WINCE Flight #97-45, 
Track 6, Center Time 20:18:05 UTC. The original threshold detects some 
of the clouds. The NPS threshold increases the area of cloud detection and 
they are better defined 
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Figure 4.23. Original cloud mask (a) and the NPS cloud mask (b) are applied to MAS 


WINCE Flight #97-45, Track 6, Center Time 20:18:05 UTC. Mask 
classifications are described in the legend. 
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V. MODIS DATA 


Figure 5.1 shows Terra’s eastern Arctic ground track for April 1, 2000. Granule 
17:00 UTC is a descending, daytime, track and is used to test the two new scene types on 


MODIS data. 


3:35 





j \ 


Figure 5.1. Terra’s eastern Arctic ground tacks for April i. 2000 (SSEC, 2000): 


A. CALIBRATION ISSUES 


There are calibration problems with some of the bands on the MODIS instrument 
and are believed to be correctable. Calibration problems with only the bands that affect 
the two new scene types are discussed below. 

1. Band 26 (1.38um) 

The objective of Band 26 (1.38um) is to detect high cirrus clouds. Unfortunately, 
there is electronic cross talk with Band 26 and bands that see the surface. This allows 


surface emittance to contaminate the reflectance value for Band 26. 
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2. Band 31 (11m) 

Band 31 (11um) has optical cross talk contamination and is on the same focal 
plane as Bands 32 - 36. After 11p1m energy passes through the filter, some of the energy 
does not reach the detector but remains in the substrate and contaminates the remaining 
Bands. 

3. Band 35 (13.94.m) 

Band 35 (13.9um) has a scan mirror angle problem. Figure 5.2 is a graphic 
illustrating the scan mirror contamination. As the scan moves from left to right, there is 
an unexpected increase in the scan mirror angle. This leads to emittance from the scan 
mirror that increases the radiance value for Band 35. Therefore, an error is induced in 
Tp(13.9um) as the scan mirror tracks from left to right from a correct value to a warmer 


value. 


Scan Track 









Correct Too Warm 


\ Granule 





Along Track 
Figure 5.2. Band 35 (13.9um) Scan Mirror Error. As the scan mirror tracks left to 
right, the mirror angle increases, which increases the brightness 
temperature across the scan. 
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B. HUDSON BAY 


Figure 5.3 shows the April 1, 2000, 17:00 UTC, MODIS granule of the Hudson 
Bay. Figure 5.3(a) shows Band 2 (0.865um) and Figure 5.3(b) shows Band 31 (11[1m). 
Hudson Bay is mainly cloud free with first year sea ice dominating the entire bay. Mid to 
low-level clouds are present along the west coast of the bay and extending northeastward 
between Coats and Mansel Islands. Finally, high clouds cover the Foxe Basin to the 


north and James Bay to the south. 


Figure 5.4 shows the Tp(11pum) test for the April 1, 2000, 17:00, UTC MODIS 
granule with a land mask applied to highlight the water areas. Figure 5.4(a) shows the 
original threshold, while Figure 5.4(b) shows the NPS threshold. As described in Chapter 
JI, the original MODIS Tp(11pm) test threshold is too high for a polar ocean region. 
When the NPS threshold is applied to the MODIS Tx,(11um) test, the clear confidence 


level for an individual pixel became more representative for a polar ocean scene. 


At the location indicated by A in Figure 5.4, the NPS threshold improves the 
detection of the clear skies in the Davis Strait because the new threshold takes into 
account the lower surface temperatures associated with polar oceans. At B in Figure 5.4, 
the NPS threshold raises the clear confidence level over the sea ice in the Hudson Strait 
where some low and high clouds are present. Finally, at location C in Figure 5.4 the NPS 
threshold raises the clear confidence level over the sea ice where the low clouds are 
present. The increase in the clear confidence level at locations B and C is not enough to 


incorrectly identify the locations as clear. 
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Figure 5.5 shows the Tp(3.9um) - Tg(11um) test for April 1, 2000, 17:00 UTC, 
MODIS granule. Again, Figure 5.5(a) shows the original threshold and Figure 5.5(b) 
shows the NPS threshold. As with the MODIS Tx,(11um) test, the MODIS T,(3.9um) - 
Tg(11um) test threshold is too high for a polar ocean region. Locations A and D are 
regions where high clouds cover the surface. The NPS threshold lowered the clear 
confidence level for these regions enough for the high clouds to be classified correctly. 
For low clouds, locations B and C had the greatest improvements. The new threshold 
detected the low level cloud that was missed at B and increased the detection of low 


clouds at C. 


Figure 5.6(a) shows the original cloud mask. It classifies portions of the clear 
area in the Davis Strait (A) as Undecided to Confident Clear. Figure 5.6(b) shows the 
NPS cloud mask correctly classifies location A as Confident Clear. This is a direct result 
of lowering the Tp(11pm) test threshold. The NPS cloud mask also changed the 
classification of the clouds at location B so that almost all of the pixels are classified as 
Cloudy. At location C, where high clouds covered James Bay, the NPS cloud mask 
correctly increased the number of pixels classified as Cloudy. The changes at B and C 
were a direct result of lowering the Tp(3.9um) - Tg(1lum) test threshold. The 
demarcation between Undecided and Confident Clear in the center of both cloud masks 


(under the label “Foxe Basin” and “Hudson Bay”) is due to Band 35 scan mirror error. 
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CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


Initial optimization of the MODIS cloud mask for the polar ocean regions is 
described based on 8 daytime MAS flights. 18 tracks (approximately 350km x 36km) 
were selected from 3 of these flights that had the greatest scene variability. MAS flights 
revealed that daytime thresholds for Tg(1 1pm) and Tp(3.9um) - Tp(11um) tests are too 
high for the polar ocean regions and must be lowered to compensated for lower average 
surface temperatures. In addition, the Tg(1lpzm) - Tp(12pm) test gives little useful 
information during the daytime. Finally, Tp(13.911m) test was never a factor in the cases 


examined in this study. 


Because of the unique nature of the polar ocean regions, two new scene types 
were developed for the MODIS cloud mask. The new scene types incorporate new 
thresholds appropriate for polar ocean regions and are included in the NPS modified 
MAS/MODIS cloud mask (NPS cloud mask). The first new scene type is Polar Ice Day 
and is initiated during the day and when the pixel is identified as having sea ice. The 
second scene type is Polar Ocean Day and is initiated during the day and when the pixel 


is identified as water located greater than 60° north and south. 


The modifications greatly improve the detection of clouds over cold polar oceans 
where sub-pixel sea ice may be present or water temperatures might falsely indicate 
clouds. The number of Cloudy pixels (< 0.66 clear confidence level) for a given scene 


was increased 12.1% on average for MAS cases. The NPS cloud mask also classified 


Eo 


two times more Probably Clear and Undecided pixels than the onginal mask due to 


greater sensitivity to thin, small clouds. 


The improvements that were observed in the MAS data were confirmed when the 
NPS cloud mask was applied to 5 uncorrected MODIS cases. One case is presented here. 
The NPS cloud mask improved the clear confidence level when surface temperatures 
were below 273°K. It also lowered the clear confidence level when there were high and 


low clouds over the sea ice. 


The NPS MODIS cloud mask that was derived from Ackerman et al. (1997) will 
provide the NIC with the means to assist its analysts by giving them a tool to correctly 
identify clouds that might have been misidentified as sea ice in the past. A better quality 
sea ice analysis will result. Ultimately, this cloud mask will also lay the foundation for 
the development of ice recognition algorithms. 


B. RECOMMENDATIONS 


The operational MODIS cloud mask must have new scene types that are 
optimized for the polar ocean regions. It is not appropriate for thresholds, set for the 
globe, to be applied to the polar oceans due to their distinctly different surface 
temperatures. Unfortunately, further work should not proceed until calibration is 
completed for the MODIS instrument. Once completed, the two new scene types must be 
tested on a larger daytime dataset consisting of Arctic and Antarctic granules. Finally, 
work needs to be done on thresholds for nighttime cases that were not possible for this 


thesis. 
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